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We report on a new search for two-body baryonic decays of the B meson. Improved sensitivity 
compared to previous Belle results is obtained from 414 fb _1 of data that corresponds to 449 million 
BB pairs, which were taken on the T(4S) resonance and collected with the Belle detector at the 
KEKB e + e _ collider. No significant signals are observed and we set the 90% confidence level upper 
limits: B(B° -> pp) < 1.1 x 10~ 7 , B(B° -> AA) < 3.2 x 10" 7 and B(B + -> pA) < 3.2 x 10" 7 . 
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Recent observations of B meson decays into several 
charmless three-body baryonic final states show branch- 
ing fractions around 10 -6 and a peak in the baryon- 
antibaryon mass spectra near threshold jl| ■ Further in- 
vestigations of the angular correlations for events in the 
threshold region Q favor a fragmentation model [H • In 
contrast, charmless two-body baryonic B decays are ex- 
pected to have lower branching fractions. However, it is 
challenging to perform conclusive theoretical calculations 
of baryon formation in B decay. Previous searches 0, H| 
for two-body decays yielded upper limits on the branch- 
ing fractions of (3 — 7) x 1CP 7 , which are consistent with a 
recent calculation using the pole model 6] but in contra- 
diction with a calculation based on QCD sum rules Q- 
Moreover, the upper limit for B° — > pp is also consistent 
with simple scaling of the measured branching fraction 
for B° — > A~p [8] by the current estimate of \V u b/V c b\ 2 , 
which gives B(B° -> pp) ~ 2.0 x 1CT 7 . 

Two-body baryonic decays provide valuable guidance 
to improve the understanding of quark/gluon fragmen- 
tation of the B meson to baryons. In this article we 
report searches with better sensitivity for the charmless 
two-body baryonic B decays B° — > pp, B° — > AA, and 
B + — ► pA The analysis is based on a 414 fb _1 data 
sample, corresponding to 449 xlO 6 BB pairs, accumu- 
lated at the T(4S) resonance with the Belle detector at 
the KEKB [l(| asymmetric e + e~ collider. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50- layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), time-of- flight 
scintillation counters (TOF), and an array of CsI(Tl) 
crystals, all located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect K\ 
mesons and to identify muons. The detector is described 
in detail elsewhere [111 . Two different inner detector con- 
figurations were used. For the first sample of 152 million 
BB pairs, a 2.0 cm radius beampipe and a 3-layer sili- 
con vertex detector (SVD1) were used; for the latter 297 
million BB pairs, a 1.5 cm radius beampipe, a 4-layer 
silicon detector (SVD2) .12] and a small-cell inner drift 



chamber were used. 

Primary charged tracks associated with candidate B 
decays are required to satisfy the following criteria: the 
track impact parameters relative to the run-by-run inter- 
action point (IP) are required to be within ±2 cm along 
the z axis (oriented antiparallel to the positron beam) 
and within ±0.05 cm in the transverse (xy) plane. 

Proton candidates are selected based on the likelihood 
functions L p , Lk and for protons, kaons, and pions, 
respectively, which are determined from particle identi- 
fication information from the CDC (dE/dx specific ion- 
ization), the ACC (Cherenkov radiation pulse height), 
and the TOF (time relative to the beam bunch crossing). 
Charged tracks with L p /(L p + L^) > 0.6 are identified 
as protons and L w /(L p + L^) > 0.6 as pions. For proton 
candidates that originate directly from B decays, an ad- 
ditional requirement, L p /(L p + Lk) > 0.6, is applied to 
improve the purity. The proton identification efficiency 
with the tighter requirements is 77% for 2 GeV/c pro- 
tons. 

A candidates are reconstructed from oppositely 
charged pion-proton pairs, satisfying the following crite- 
ria: the separation distance between the pion and proton 
at the A decay vertex must be less than 12.9 cm along 
the z axis; the distance of the closest approach in the xy 
plane to the IP is greater than 0.008 cm for each track; 
the flight length of the A candidates must be greater than 
0.22 cm in the xy plane; the angular difference between 
the A momentum vector and the vector between the IP 
and the A decay vertex must be less than 0.09 rad. Fi- 
nally, the reconstructed invariant mass of A candidate is 
required to be within the mass interval (1.116 ± 0.005) 
GeV/c 2 , corresponding to ±3<r in the mass resolution. 

B candidates are reconstructed by pairing two 
baryons and are identified with two kinematic vari- 
ables: the beam-energy constrained mass, Mb c = 
V -^bcam/ c4 — Pb I ° 2 i an d the energy difference, AE — 
Eb — -Ebcam, where i^bcam is the run-dependent beam en- 
ergy in the center-of-mass (cm.) frame and pb and Eb 
are reconstructed momentum and energy of the recon- 
structed B candidates in the cm. frame, respectively. 
The signal region is defined as 5.27 GeV/c 2 < Mbc < 
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5.29 GeV/c 2 and \AE\ < 0.05 GeV, while the sideband 
region is defined as 5.2 GeV/c 2 < M hc < 5.26 GeV/c 2 
and \AE\ < 0.2 GeV. 

All possible backgrounds are investigated using a 
GEANT3 [H based Monte Carlo (MC) simulation. 
Backgrounds from the charmful and charmless B decays 
are found to be negligible. The dominant background is 
from continuum e + e~ — ^ (?Q (<7 = u, d, s, c) events, which 
are studied using the sideband data. 

Continuum background is suppressed by requiring 
| cos#th| < 0-9, where t h is the angle in the cm. frame 
between the direction of one B daughter and the thrust 
axis 1J| , formed by the particles not associated with the 
B candidate. The cos#th distribution is nearly flat for 
signal but strongly peaks at ±1 for the continuum. 

The background rejection is further improved using 
event topology, B candidate vertex and B flavor tagging 
information. First, we combine a set of modified Fox- 
Wolfram moments [l5| into a Fisher discriminant [l6| 
T to distinguish spherically distributed BB events from 
the jet-like continuum backgrounds. Figure [Ha) shows 
the Fisher discriminant for pp signal MC and sideband 
events. We then use the cosine of the angle 9b between 
the B candidate flight direction and the z axis (Fig.QJb)). 
For the pp mode only, we use the distance Az along the 
z axis between the pp vertex and the vertex formed from 
the remaining charged tracks (Fig. [He)). The Az prob- 
ability density function is modeled by a double Gaus- 
sian independently for the SVD1 and SVD2 data due to 
the better vertex resolution provided by the SVD2. The 
PDFs for cos 9 b and T are described by a 2nd order poly- 
nomial and a double asymmetric Gaussian, respectively. 
The quantities T , cos 9b and Az (for the pp mode only) 
are combined to form likelihoods Cs and Cb for signal 
and continuum background, respectively. The normal- 
ized ratio TZ = Cs/(Cs + Cb), shown in Fig. QJd), peaks 
at unity for signal and at zero for continuum. 

The distribution of TZ is somewhat correlated with the 
event's B flavor information. The standard Belle flavor 
tagging algorithm [ljj provides a tagging quality factor 
r that ranges from zero for no flavor tagging information 
to unity for unambiguous flavor assignment. We apply a 
requirement on the likelihood ratio TZ depending on the 
value of r (and on the inner detector configuration for 
the pp mode) . The likelihood ratio requirements listed in 
Table [I] are determined by optimizing the figure-of- merit, 
Ns/y/Ns + Nb, where Ns and Nb are the expected sig- 
nal and background yields. The signal yields are esti- 
mated from the product of the number of BB events, 
the signal efficiency from MC and the assumed branch- 
ing fraction of 10 -7 . The expected background yields 
are obtained by scaling the amount of data in the region, 
\AE\ < 0.05 GeV and 5.2GeV/c 2 < M hc < 5.26GeV/c 2 , 
to the signal region. The Mb c distribution of the con- 
tinuum is assumed to be independent of AE. Therefore, 
we model the M^ c distribution with an ARGUS func- 
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FIG. 1: Distributions of (a) the Fisher discriminant, (b) 
cos 9b, (c) Az, and (d) likelihood ratio for B° — » pp can- 
didates . The solid histogram is for signal MC, while the 
dashed histogram is for sideband events. 



tion using the data in 0.1 GeV < \AE\ < 0.2 GeV. 
The scaling factor is estimated to be the ratio of the area 
of the Mbc signal region to the sideband region based on 
the ARGUS function. 



TABLE I: Minimum values of the likelihood ratio TZ for three 
ranges of the flavor tagging quality r (and for the inner de- 
tector configuration for the pp mode). 





< r < 0.5 


0.5 < r < 0.75 


0.75 < r < 1.0 




SVD1 SVD2 


SVD1 SVD2 


SVD1 SVD2 


B° -> pp 


0.85 0.85 


0.8 0.65 


0.65 0.65 


B° -> AA 


0.85 


0.6 


0.35 


B+ -> P A 


0.8 


0.75 


0.65 



The B signal efficiency for each mode is obtained us- 
ing the MC simulation after applying all analysis require- 
ments except the proton identification. The B signal effi- 
ciency is corrected for the proton identification; the pro- 
ton identification efficiency is estimated using A — > p7r~ 
decays in data. The systematic error in the B signal 
efficiency arises from proton identification, tracking effi- 
ciency, A selection and the likelihood ratio requirement 
in each flavor-tagged region. The statistical error in 
the proton efficiency obtained in the A sample is in- 
cluded in the systematic error for proton identification. 
The tracking systematic error is studied using an inclu- 
sive D* + — > D°tt + sample, where D° — > Kgir + ir~ and 
Kg — > 7r + 7r~. Candidate D* + mesons are kinematically 
identified using the momentum of the slow pion from the 
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D* + decay, D° and Kg mass constraints, and the tra- 
jectory of one of the pions in the K s decay. The track- 
ing efficiency is measured by searching for the other K s 
daughter; by comparing the efficiency in data with the 
Monte Carlo expectation, the tracking systematic error is 
found to be 1% per track for p > 0.3 GeV/c but slightly 
larger for low momentum particles. The systematic error 
associated with requiring a detached vertex for the A can- 
didate is studied by comparing the ratio of D + — > K s n + 
and D + — ► K~ir + Tr + yields with the MC expectation. 
The resulting Kg detection systematic error is ±4.5%, 
which is used as the A systematic error. The likelihood 
ratio requirement is studied using control samples with 
the same number of charged particles in the final state: 
B° -> K+ir- for pp, B° -> D~>k + -> (If +7r-7r-)7r+ for 
AA and B+ -> lfn + -> (K+tt-)tt+ for pA. The total 
systematic uncertainty is computed by adding the cor- 
related errors linearly, and then adding the uncorrelated 
ones in quadrature. The systematic errors are summa- 
rized in Table El 



TABLE II: Summary of systematic errors, given in %. 





B° -> pp 


B° -> AA 


B+ -> pA 


Tracking 


2.00 


4.29 


3.16 


PID 


0.40 


0.12 


0.26 


72. Requirement 


3.70 


1.22 


0.75 


A Selection 




4.50x 2 


4.50 


B(A -> pK~) 




1.56 


0.78 


# of BB 


1.27 


1.27 


1.27 


Total (%) 


4.42 


10.24 


5.76 



The expected background contribution is obtained us- 
ing the sideband region described above. A loose 1Z re- 
quirement is applied to ensure large statistics to deter- 
mine the ARGUS parameters. The systematic uncer- 
tainty in the backround prediction comes from the fit er- 
rors of the ARGUS parameters and the statistical errors 
of the events in the sideband region. 

The estimated background yields in the signal region 
are listed in Table Ifffl and are close to the number of 
observed events. Since no statistically significant signals 
are found, we follow the Feldman-Cousins approach [l9| 
to estimate 90% confidence level (C.L.) upper limits on 
the signal yields, using the implementation of J. Conrad 
et al. 20] to include the systematic errors. The final 
results are listed in Table Ifffl 

In summary, we have performed a search for the decays 
B° — > pp, AA, and B + — > p A in a sample of 449 x 10 6 
BB events, which is three times larger than the dataset 
used in the previous analysis 0] ■ We find no evidence for 
signals and place 90% C.L. upper limits on the branch- 
ing fractions of 1.1 x HT 7 , 3.2 x 10~ 7 , and 3.2 x 10~ 7 
for the pp, AA, and pA modes, respectively. These up- 
per limits improve our previous results 0] and are more 
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FIG. 2: _Scatter plots of AE vs._M b c for B° -> pp (upper), 
B° —t AA (middle), and B + — > pA (lower). The signal region 
is indicated by the rectangle in each plot. 



TABLE III: Summary of the B° -> pp, AA and B+ -» pA 
searches, where e is the reconstruction efficiency including the 
sub-decay branching fraction, iV bs is the observed number of 
events in the signal region, N^ p is the expected background in 
the signal region, Nqo is the yield limit at the 90% confidence 
level, and BF is corresponding upper limit for the branching 
fraction. The uncertainty in is the systematic uncer- 

tainty due to the ARGUS parameters and statistical error of 
the sideband sample. 



Mode 




e [%] 


N obs 


N bg 


^90 


BF [icr 7 ] 


B° - 


PP 


17.73 


25 


26.8±2.3 


8.7 


< 1.1 


B° -> 


AA 


4.24 


3 


1.2±0.5 


6.1 


< 3.2 


B+ - 


• pA 


8.44 


16 


12.9±1.7 


12.1 


< 3.2 



stringent than other experimental limits [5j, l2l| . More- 
over, although our pp and pA results already reach the 
pole model predictions 0], no clear signals are seen. 
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